A Stochastic Description for Extremal Dynamics 
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We show that extremal dynamics is very well modelled by the "Linear Fractional Stable Motion" 
(LFSM), a stochastic process entirely defined by two exponents that take into account spatio- 
temporal correlations in the distribution of active sites. We demonstrate this numerically and 
analytically using well-known properties of the LFSM. Further, we use this correspondence to write 
an exact expression for an n-point correlation function as well as an equation of fractional order for 
interface growth in extremal dynamics. 



minimum pinning force (and hence closest to its instabil- 
ity threshold) advances first, thus h{xo) h{xo) + Ah. 
At this new position, a new random pinning force is en- 
countered /p(xo,t -I- 5t) = fo{xo, h{xo) + Ah). The exter- 
nal loading F on the system, and interactions along the 
front, produce a local driving force on each site x propor- 
tional to f{x, t) = F. J G{x—y)h{y, t) dy where the kernel 
G{x) (X |a;|~''~^ accounts for long range interactions me- 
diated by the medium. The loading F is then adjusted 
so that only one site depins f{x,t) — fp(x,t)] the others 
remain trapped since f{y,t) < fp{y,t) for y ^ x. The 
dynamics of advancing the minimum site and readjusting 
the others is continued indefinitely. 

A wide class of extremal models have already been 
studied extensively by Paczuski et al. 0. These mod- 
els include the Bak-Sneppen evolution model |^ and the 
Sneppen Interface Model [0. All these models try to 
explain driven motion under strong pinning by means of 
a discrete, deterministic dynamics. Only one site is ac- 
tive at every instant of time. However the "time" is only 
a way to index the sequence of events. Further the ex- 
tremal condition can be thought of as a way to retain 
the information of the spatially quenched heterogeneities 
that determine the evolution of the front. All the infor- 
mation in this sort of dynamics is clearly contained in 
the "activity" map: a space-time plot of where the front 
is active at every instant of time. 

Previous studies regarding extremal models P,pd| have 
shown that most of the relevant information is contained 
in the probability density function (hereafter, pdf) of the 
activity map. In the stationary regime, assuming that 
the activity was located at xq at time the probability 
that it is at x at time t is: 



Since the eighties, several models have been proposed 
10-^ to describe in a common language diverse physical 
situations such as roughening of a crack front in frac- 
ture jl], wetting front motion on heterogeneous sufaces 
1^, dynamics of a ferromagnetic domain wall driven by 
an external magnetic field Q , motion of vortices in type- 
II superconductors [Q , fluid invasion in porous media , 
solid friction |g] or even biological evolution 0. All these 
models propose to explain the dissipative behaviour of 
the system by the competition between an elastic restor- 
ing force and a non linear, randomly distributed, time- 
independent, pinning force. In the case of the spreading 
of a partially wetting liquid for example, the pinning force 
is due to surface chemical heterogeneities or roughness, 
and the elastic restoring force is a result of the surface 
tension at the liquid/vapor interface. In the case of strong 
pinning it is well known |^ that, when subjected to an 
external driving, the wetting front displays local instabili- 
ties that force it to advance quasistatically. This property 
makes it difficult to handle the problem in the continuum. 
In fact, in the stationary regime, the main contribution 
to the global displacement is from jumps of local parts of 
the chain resulting from these instabilities Recently, 
Tanguy et al. Q have proposed to describe this sort of 
evolution by an extremal model: only the site closest to 
its instability threshold advances. After a jump, the in- 
stability thresholds of all the sites are modified by the 
(elastic) couplings between sites. More precisely, in their 
model (hereafter LREM for Long Range Elastic Model), 
the interface of size L is defined on a discrete lattice 
{x, h). Initially the front h{x) = 0, and the pinning forces 
fp{x,t = 0) = fo{x,h = 0) are assigned independently 
from a flat distribution. The site xq subjected to the 
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i( \ ( r""~i for r > 1 /r,N 
|_ r' for r ^ f . 

While the exponent a controls the asymptotic behavior 
of the time independant function (jj, z controls the prop- 
agation of the activity along the system as a function 
of time and is therefore termed the "dynamical expo- 
nent" . The above distribution is self-afhne and therefore 
its temporal evolution is completely defined through the 
exponents z and a Moreover, as can be seen from 
Eq. the form of the distribution of \x{t) — xo\ is in- 
dependant of t for large enough t: it is "a-stable" [^s). 
It has been shown in that a — b. It is hence also 
possible to consider the activity as performing a "Brow- 
nian like" motion similar in spirit to studies of anomalous 
diffusion The relevant parameters here are: the ex- 
ponent a appearing in the stationary distribution of the 
distance between successive active sites (distribution of 
increments) : 

p{{x{t + l)-x{t)) = l)(xl-''-\yt (3) 

and the exponent H , characteristic of the moments: 

{\xit) - x{t')\''y^'' w |t - t'\" for a<a. 

H accounts for possible temporal statistical dependence 
between jumps. When there is no temporal correlation 
between jumps, this distribution is Brownian for a > 2 
and a "Levy flight" for < a < 2. In the former case 
H = 1/2 and in the latter case it is easy to show, for 
example from the asymptotic expression of p{x,t) [ p^ , 
that H — l/a. However when there are temporal corre- 
lations, as in some of the models studied in or in the 
case when elastic interactions are long ranged 1^, this is 
no longer true and H ^ l/a. The value of H is hence 
indicative of the presence or absence of temporal corre- 
lations in the jumps of the activity. Attempts have been 
made in the past to understand the space-time plot of 
the activity in extremal models as an uncorrelated Levy 
flight 1^,^ , by keeping only the exponent a and thus as- 
suming H — l/a. However, numerically, this hypothesis 
leads to erroneous quantitative predictions. Thus there 
are long range time correlations which have to be incor- 
porated in the description. In this Letter, we introduce 
and study such a model. 

We define here a Linear Fractional Stable Motion. Let x 
denote a process generated in the following manner: 

x(t)=^/(t,u)r;(«), (4) 

where ri{u) is an uncorrelated noise with a symmetric 
distribution p{rj = x) ~ |a;|~"~"'^. Since we would like 
to consider stationary processes (if our eventual aim is to 
describe steady states) we assume that f{t, u) = f{t — u). 
The definition above implies basically that x consists of 
a sum of uncorrelated levy jumps weighted in time by 



f{t — u). This weight function therefore controls the time 
dependence of the statistical properties of x. It is easy 
to show that the sum in Eq. ^ can be performed much 
as for independant levy flights and the random variable 
X = x{t) (given that X = at t = 0) has a probability 
density function 

p{X,t)^ y"exp(ifcX) exp(-<|fc|")dfc/27r. (5) 

where (ctj)" = ELo ^^r the so-called LFSM 

|l3| , the function f{t,u) reads: 

f{t,u) = {t-uf 

where the parameter d satisfies —l/a<d<l — l/a. 
The LFSM is thus a self-similar process with stationary 
increments. Further the exponent H defined earlier is 
here = d + l/a and = af\t\°'^. Thus x{t) can be 
seen as an integral or derivative of fractional order, of 
the noise rj. 

In order to have a prescription for constructing the 
LFSM given the activity plot of the extremal models 
that we studied, it is important to determine carefully 
the H exponent in these models. We have done this 
by performing first a wavelet transform Jl6t , which con- 
sists of computing wavelet coefficients dx{a, k) = {x, ipa.k) 
where 4'a.k{t) = -\/aV'o[a(i — k)] is a collection of dilated 
and translated templates of the mother wavelet tpoit)- 
The wavelet transform is a relevant tool |jl^ to anal- 
yse self-similarity because the dx{a,k) of a self-similar 
process with stationary increments i) reproduce exactly 
the self-similarity, ii) form stationary sequences, iii) are 
quasi independent statistically one from another. For the 
LFSM, it has been shown |l8| that H can be estimated 
by performing a linear fit in a (log2(a),ya) plot, where 
Ya — (log2 \dx{a, k)\)k. Confidence interval for the esti- 
mate of H can be theoretically derived ||l^. Figure |l| 
illustrates this estimation procedure on data produced 
by an LREM using 6 = 1 g]. Only one such example is 
shown, but many other trials were performed using other 
values of &, or other models such as the Bak-Sneppen 
model. All the models produced power-laws of similar 
quality, thus proving that the activity map of these mod- 
els can be modelled as an LFSM. 

We now try to compare the exponents a and z pre- 
dicted by the scaling form (|l|) obtained for an extremal 
model (such as the LREM) and its LFSM description. 
From the definition of the LFSM given by Eq. (|^) and 
using the fact that erf — ai\t\"^ , it is easy to show that 
the probability density function for the LFSM satisfies 
the scaling form (|^) with z — 1/H. This is also very well 
verified on numerical simulations of LFSMs (fig.|[b) and 
gives an estimate of z with a precision of ±10^^ for a 
signal of 10^ time steps. To compare, this analysis has 
also been applied on the LREM (see fig.H.a) allowing es- 
timates on z with a precision of ±0.1. The two estimates 
of z give us the same result. 
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To further justify the fact that the LFSM models ac- 
curately the activity map of extremal models, we stud- 
ied another property, namely the distribution of return 
times, which had been characterized on LREM as power- 
laws, but never on LFSM. The first return time is the 
time T elapsed between two subsequent activities at a 
given site. It has been shown for LREM, that T is dis- 
tributed as a power-law p{T) oc T^'^f up to a maximum 
time T* such that the activity has spead over the entire 
system, T* (x . For the LFSM, we have already seen 
that the activity plot generated has a dynamic exponent 
z = 1/H. When H < 1, z is also the fractal dimension 
of the activity plot. From this it is easy to show that the 
first return time exponent (which is simply a one dimen- 
sional cut of a fractal set) should be 



Tf = 2- H = 2-l/z. 



(6) 



This is a relation already known to numerically hold 
for extremal models 0. Hence we see that the two are 
again identical in this regard. We have also checked this 
numerically for different extremal models and their cor- 
responding LFSM's. It is important to note that, when 
H > 1, TpiRST = 1 for all H . Thus our identification of 
the LFSM nature of the activity is indeed justified, and 
allows us to access a novel property of the LFSM which 
has never been reported so far in the literature. 

We now turn to some of the interesting consequences 
of having an analytical expression (such as Eq. ^ for the 
space-time plot of the activity. Consider the two-point 
function Pijijh)- the probability of having a jump li 
and I2 at two consecutive instants. Using Eq. we find 
that the expression of this function in Fourier space is just 
~ / exp{ik{h+l2)) exp(-|2|"^|fc|")dfc. In real space this 
is just the function l/ih + h)'^'^^ asymptotically. We can 
use this to compute the conditional probability P{l2/li) 
— the probability of having a jump of length I2 at time 
t — 2 given that there was a jump li at time t — 1 — 
using the well known relation P{l2/li) = Pih^h)/ P[W)- 
We see that this is just ~ 1/(1 + h/WT^^ thus the- 
oretically confirming the numerical scaling obtained in 
We can now generalise to the full n-point function 
P{h^h,h----ln) (discussed in |l^) of having a jump of 
length li at time f = 1, a jump of length I2 at time 
t = 2 and so on till t = n. Systematic expansions can be 
obtained for the conditional probabilities just as for the 
two-point function and can also be verified numerically 
on the models. The theoretical understanding of these 
correlation functions is hence one of the achievements of 
our mapping. 

A further very interesting consequence of the definition 
of the LFSM is the following equation that we can write 
down for the front propagation under this dynamics. To 
do this, we define the height of an interface at time t at 
a spatial location X as simply the accumulated activity 
there upto time t: 



h{X,t)=^Ah.[ 5{X ~ x{t'))dt' . 



(7) 



Note that we can use with this definition (0) to perform 
the scale transformation commonly used for self-affine 
surfaces: X bX , t ¥'t and h — * b^h, where x is 
the so called roughness exponent for the height. A power 
counting on both sides of this equation gives the relation 
X — z—\, known to hold for extremal dynamics . 

Using the definition of the LFSM (jj) and the following 
definition of the fractional derivative lll9|: 



{X^x{t))^ 



1 d r^'^+ U{u)-^{x{t)) 

r(i - /3) i™o ds {x{t) + 5~u)P 

we obtain the following equation for the height after te- 
dious calculations (the details of which we leave for a 
longer paper [po|): 



d""h{X,t) _ T{aH + l) 



T{2-aH) 



du.rjiu). 



r(a + l) 

'^^ dXdu' 



(9) 



In deriving the above, we have basically used some of the 
techniques developed for derivatives of fractional order 



such as the fractional Taylor expansion 1 19 1 and fractional 
integration by parts, c is a constant with c ^ {51)^^°"^ 
taking into account the right dimensionality of the noise. 

Wc have thus changed a problem with quenched disor- 
der to one with a multiplicative annealed noise, rj{t), with 
our approach. There are several points worth mentioning 
in this regard. Using the fact that p{X, t) is proportional 
to the number of times the site X has been visited at 
time t (for different realisations of the process), that is 
p{X,t) ~ {dh{X,t)/dt), we can also conclude from Eq. 
(^ that an effective Fokker Planck equation for the ac- 
tivity should be one of fractional order in space and time: 



d'^"piX,t) ^.y,. 



(10) 



We note that the expression Eq. (^ is consistent with 
the above equation in the long-time limit [pT| . For an 
uncorrelated process the above equation reduces to the 
one studied in |2^. An equation such as the above has 
also been derived for a process with long-ranged jumps 
and a power-law waiting time between jumps p^ ]. It 
would be interesting to see the connection between the 
above process and Eq. (^) more microscopically. 

There are also several other interesting points to in- 
vestigate. As we have proved above, any extremal model 
can be modelled by a tangent process, which we call the 
LFSM, for a suitable value of H and a. In Fig. ^ we have 
indicated the above parameter values for all the extremal 
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models mentioned in this letter. It is interesting to note 
that for the models studied in ||] a further relation ex- 
ists between H and a; H = 3/2{a +1) 1^ as indicated 
in Fig. 3. The two other extremal models studied (The 
Bak Sneppen model and the Sneppen Interface model) 
however do not seem to obey this relation (Fig. 3). It 
would be interesting to understand why the above rela- 
tion exists for the LREM. 

In summary we have proposed the LFSM as a stochas- 
tic model for extremal processes. We argue that the self- 
affinity and a-stability of this process, together with the 
stationarity of its increments makes it an accurate de- 
scription of extremal dynamics. We demonstrate this by 
first determining the value H and a for all the extremal 
models listed above by using a wavelet analysis of the ac- 
tivity. Subsequently we show that a LFSM process gener- 
ated using these values shows the same scaling predicted 
by Eq. |l| because of a simple relation between the dynam- 
ical exponent z and the exponent H. In order to test such 
a correspondence further we also investigate a property 
which has never been studied in the LFSM, namely the 
time interval distribution between recurrences of activity 
at a particular site. Further the LFSM is amenable to 
an analytical treatment much more easily than extremal 
dynamics. We show how the scaling form (|l|) is a simple 
consequence of the definition of the process. We comment 
on how the process allows us to understand the full n- 
point probability distribution for the increments. Finally 
we conclude by proposing a fractional partial differential 
equation for a front subjected to extremal dynamics. 
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FIG. 1. Log-log plot of the amplitude of the wavelet com- 
ponent versus scale factor using a DaubechiesS wavelet on 
time series produced by a LREM [6]. It exhibits a quasi per- 
fect power- law revealing self-afBnity. The slope is H + 1/2 
where H is the self-affinity parameter. 
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FIG. 2. Data collapse of eight different probability distri- 
butions p{x, At) with time intervals At ranging from 1 to 128 
(a) for the extremal model of ^(L=16384). a — 1.45 and 
the best data collapse is obtained for z = 1.55. Idem for a 
LFSM (b) with a = 1.5 and H = 0.562. The best collapse is 
obtained with z — 1/H. 
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FIG. 3. This diagram represents the (wavelet-based) esti- 
mated Hs (and their confidence intervals) for various a and 
various extremal models. The dot-dash curve 1/a is the esti- 
mate for an uncorrelated levy-flight. The function 3/2(a -I- 1) 
is the numerically found best estimate for z in the LREM. 
The two other extremal models studied however do not seem 
to obey this relation. 
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